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^ : ABSTRACT 

1^ ■ We report on the feasibility study of a Laue lens for hard X-rays (> 60 keV) based on mosaic crystals, for 
astrophysical applications. In particular we discuss the scientific motivations, its functioning principle, the 
procedure followed to select the suitable crystal materials, the criteria adopted to establish crystal dimensions 
I and their distribution on the lens in order to obtain the best lens focusing capabilities, and the criteria for 
^SJ ' optimizing the lens effective area in a given passband. We also discuss the effects of misalignments of the crystal 
K>" ' tiles due to unavoidable mechanical errors in assembling the lens. A software was developed to face all these 
topics and to evaluate the expected lens performance. 
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'st : 1. SCIENTIFIC MOTIVATIONS 

o 

The role of hard X-ray astronomy (> 10 kcV) is now widely recognized. The numerous results obtained with 
the most recent satellite missions {BeppoSAX , Rossi-XTE) on many classes of X-ray celestial sources have 
demonstrated the importance of the broad band (0.1 300 keV) spectroscopy in order to derive an unbiased 
O \ picture of the celestial source physics, like to establish the source geometry, the physical phenomena occurring in 
■ the emission region, the radiation production mechanisms, an unbiased separation of the contribution of thermal 
^ ] emission phenomena from the phenomena due to the presence of high energy plasmas (thermal or not thermal) 
!• ■ and/or magnetic fields and/or source rotation. 

■ ""j In spite of the excellent performance of the high energy instrument PDS (Phoswich Detection System) aboard 

rN BeppoSAX ,^ the most sensitive instrument ever flown in the 15-200 keV energy band, even in the case of the 
strongest Galactic (e.g., Cyg X-1 in soft state, Her X-1) and extragalactic (e.g., 3C373, MKN 3) X-ray sources, 
the statistical quality of the measured spectra becomes poor in the highest part of the instrument passband 
(> 80 keV). Thus the development of focusing optics in this band and, more generally, in the entire passband 
covered by the BSAX/PDS and possibly beyond it, is of key importance to overcome the limitations of the direct 
viewing telescopes (with or without masks) and to allow the study of the high energy spectra of the celestial 
sources with the same detail which is achieved at lower energies (< 10 keV), where focusing optics are available. 

In fact. X-ray mirrors, based on the external reflections with very high focal lengths (> 50 m), or 'supermir- 
rors' based on Bragg diffraction from multilayers of bi-strates made of high and low Z materials (e.g., Joensen et 
al.^), can overcome the sensitivity problem up to about 70 keV. For higher energy photons, an efficient focusing 
is a much more challenging task. 

Goal of our project is the development of a focusing telescope which efficiently focus hard X-/gamma-rays in 
a broad continuous band, from 70 keV to > 300 keV, by exploiting the Bragg diffraction from mosaic crystals in 
Laue configuration. 
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2. LAUE LENS CRYSTAL REQUIREMENTS 



The X-ray diffraction from crystals in reflection (Bragg geometry) or transmission (Lauc geometry) confignrations 
is well known (e.g., James*^). The first step of our feasibility study we are reporting on is to establish the lens 
physical properties, like the most suitable crystal materials, their sizes (cross section, thickness), and their 
orientation in the lens, which maximize the lens performance with sustainable costs. 

In order to better understand the problematics of a Laue lens, we first summarize basics of crystal diffraction 

and properties of mosaic crystals. 

2.1. Basics of the crystal diffraction theory 

As it is well known, the X-ray diffraction from crystals is based on the two equivalent equations: 



the first due to M. von Laue,'* and the second due to W.H. and W.L. Bragg. ^ 

In the Laue equation (1), kj and kj are the wave vectors of the incident and the diffracted beams and g is 
the reciprocal lattice vector. In the Bragg law (2), dhki is the distance between the lattice planes with Miller 
indices hkl, 9b is the scattering angle for diffraction (the Bragg angle), n is the diffraction order and A is the 
wavelength of the photon. The energy E of the photon is related to A through the relation E = hc/X, in which 
he = 12.4 keV-A. The two equations are completely equivalent and provide the same information: for example, 
the angle between the two wave vectors ki and k<i is exactly twice the Bragg angle 6b and the reciprocal lattice 
vector g is perpendicular to the planes and its module has magnitude 27rd^^j. 

According to the above equations and depending from the crystal structure, a parallel polychromatic photon 
beam will be splitted in several monochromatic beams with directions depending on the photon energies present 
in the beam and on crystal planes which diffract photons. In fact, due to crystal defects, ion thermal motion 
and finite size of the sample, the diffracted beam is not monochromatic but shows a small spread in energy with 
an intensity peak centered around the energy which corresponds to 6 through the Bragg law. 

We have to face a couple of important problems for getting the Laue lens for our applications: a geometrical 
problem, and an efficiency problem. The first concerns the need of focusing photons into a very small area, the 
second is the need of maximizing the reflection efficiency in a broad band. A nearly-perfect crystal, could solve 
the flrst problem through a proper geometry of the crystal shape, but not the second one (the high efficiency 
is limited to a very narrow band). The best compromise is achieved by using imperfect crystals with a suitable 
spread of the photon energies which can be reflected. The best crystals which can satisfy our requirements are 
the mosaic crystals, as we demonstrate below. 

2.2. Mosaic crystals 

Mosaic crystals are made of many microscopic perfect crystals {crystallites) with their lattice planes slightly 
misaligned with each other around a mean direction, the one corresponding to a reciprocal lattice vector go. 
In the configuration assumed, the mean lattice plane is normal to the mosaic crystal surface. The distribution 
function of the crystallite misalignments from the mean direction can be approximated by a Gaussian function: 



where A is the magnitude of the angular deviation from the mean. The Full Width at Half Maximum (FWHM) of 
the Gaussian function defines the mosaic spread (3 oi a mosaic crystal. The quantity W{A)dA is the probability 
of finding crystallites with planes misaligned by an angle in the range [A, A + dA] from the mean direction. 

For the Laue geometry and plane parallel plates, with diffracting planes perpendicular to the plates (see the 
right side of Fig. 4), the intensity of the diffracted beam Id{A, E) is given by the following equation^: 



ki - kd = g 



(1) 



2dhM sin 0B =nX, 



(2) 




(3) 



Id{A, E) = lo sinh (aT) exp - (/x + 7o(t) — = (1 - e-^^'^)e-< 



(4) 



where Iq is the intensity of the incident beam, /j, is the absorption coefficient per unit of length corresponding 
to that energy, 70 is the cosine of the angle between the direction of the photons and the normal to the crystal 
surface, T is the thickness of the mosaic crystal and a is: 



a = a{E,A) = W{A)Q{E) 



(5) 



where 



Q{E) = 



reFhki ' 3 l+cos^(2^ 
V 2sm26B 



(6) 



in which re is the classical electron radius, Fhki is the structure factor, V is the volume of the crystal unit cell, 
A is the wavelength and 0^ is the Bragg angle for that particular energy. 

The quantity 7ocr is known as the secondary extinction coefficient and T/70 is the distance travelled by the 
direct beam inside the crystal. 

Prom Eq. (4) it is possible to see that, for fixed values of T, the value of Ij, is monotonically increasing with 
(T and monotonically decreasing with fi; thus to increase the diffracted beam, it is necessary to adopt materials 
with low values for /x and great values for a. The best thickness, Tb^gt, that maximizes the value of Id can be 
obtained equating the first derivative of (4) with respect to T to zero. The result of this operation is: 



The maximum achievable value for 7^ is 1/2, that is a limiting value obtainable in the ideal case of zero absorption 
and infinite secondary extinction. Using the Eq. (4), it is possible to deduce some important and useful properties 
of the diffracted photon beam. 

2.3. Properties of the photon beam diffracted by mosaic crystals 

The total diffracted intensity of a parallel photon beam increases with the mosaic spread p. However also its 
divergence increases with /?. As a consequence, a polychromatic parallel beam, after its reflection by a mosaic 
crystal not only has a different energy distribution but also a different size. We have determined the spread of a 
photon beam diffracted by a mosaic crystal in a plane perpendicular to the incident beam direction. Assuming 
a crystal with plane parallel plates and sizes shown in Fig. 4, with mean lattice plane perpendicular to the 
crystal surface and misalignment distribution of the microcrystals (see Eq. 3) along all azimuth directions, for 
a uniform parallel photon beam incident on the crystal, which is large enough to intercept the whole of it with 
incidence angle 6b (grazing angle from the mean lattice plane), the diffracted beam in a plane perpendicular to 
the incident beam has dimensions larger than the crystal cross section. The maximum enlargement is along the 
direction perpendicular to the mean lattice plane. In this direction the linear enlargement (where is 98% of the 
reflected power) is 2F x (3, where F is the distance from the crystal, while in the perpendicular direction it is 
F t;in{20B)P- Given the small Bragg angles for hard X-rays (> 60 keV), the latter quantity is a factor 100 lower 
than the former. 

As far as the energy distribution of the diffracted beam is concerned, assuming a uniform polichromatic 
beam, by means of the Bragg law (see Eq. 2), we get that the energy range diffracted by the mosaic crystal 
along the direction normal to that of the mean lattice plane is approximately given by (i?i, E2), where Ei = 
Eb{1 — (3 /{29b)), and E2 = Eb{1 + (3 /{29b)), in which Eb is the Bragg energy corresponding to the incidence 
angle 9b- Thus the energy spread corresponding to the mosaic spread /? is Eb(3/9b- Given that for hard X-rays, 
9b ~ hc/{2dhkiE), the energy spread is proportional to E"^. 

2.4. Crystal material choice 

From the formulae given in Section 2.2, the materials for a Lane lens should be chosen according to the following 
guide-lines: low ji and high Q to maximize the crystal reflectivity, low dhu to get large reflection angles, and thus 
high lens effective area (see Section 5), in the lens energy passband. However mosaic spread (3 has to be chosen 
with care. Indeed, on one side a higher /3 value gives a higher crystal reflectivity, on the other side it degrades 
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Figure 1. Energy behavior of /i and a for the materials investigated: HOPG (002), Copper (111) and Germanium (111). 
The cr values are computed assuming A = and 77 = 1 arcmin. 

the focusing capabilities of the lens. Apart from that the prerequisite for the crystal material is either that it is 
already available on the market with the required mosaic properties or that its development with the required 
mosaic properties is feasible. 

Thus we started our feasibility study assuming mosaic crystal materials that satisfy the the above prerequisite. 
These include Highly Oriented Pyrolitic Graphite (HOPG), Germanium and Copper. The lattice structure of 
HOPG (with Miller indices 002) is hexagonal with spacing (ioo2 = 3.874 A and that of Cu (111) and Ge (111) is 
Face Centered Cubic (FCC) with d\\\ = 2.084 A for Copper and dm = 3.268 A for Germanium. 

HOPG is available since a long time, but its mosaic spread is /3 > 0.3°, while Ge (111) can be obtained in 
mosaic configuration with a spread f3 < 1', and mosaic crystals of Cu (111) can be now obtained with a mosaic 
spread between 1 and 10 arcmin. Given their low mosaic spread, the last two materials show bettor focusing 
capabilities. Instead HOPG appears less suitable, not only for its larger mosaic spread, but also for an additional 
drawback for its use in transmission configuration. Indeed, given that it is produced in plates with a maximum 
thickness of a few millimeters and the (002) planes are parallel to the plates, its use in Lauc configuration with 
lattice planes normal to the crystal cross section, requires several plates to be stacked together, a task which we 
have experimented and found to introduce further mosaic spread. In spite of these technological problems, we 
have compared the expected performance of all these crystals. The energy dependence of their cr and /i in the 
60-400 keV band is shown in Fig. 1, in which for a we assumed 77 = 1 arcmin and A = 0. The energy behavior 
of the expected peak reflectivity for rj = 1 arcmin, a crystal thickness of 0.5 mm for HOPG and 0.2 mm for Cu 
(111) and Ge (111), is shown in Fig. 2. As can be seen, Cu(lll) shows the best peak reflectivity in the 100-400 
keV energy band. 

3. LENS GEOMETRY AND ITS MAIN PROPERTIES 
3.1. Lens geometry 

For a parallel photon beam, as that coming from a celestial X-ray source, among the possible focusing geometries 
of a Laue lens made of mosaic crystals in transmission configuration with mean lattice plane perpendicular to 
the crystal surface, the paraboloids of revolution or the spherical surfaces appear the simplest. They satisfy 
the projective transformations, as required for a perfect imaging. * The crystals should be curved according to 
the chosen geometry. In the case of flat crystals, the best approximation of the geometry is obtained by using 
small crystals. For high focal lengths F (> 2 m), as requested in the case of hard X-rays, the two geometries, 
paraboloid of revolution or sphere, are practically undistinguishablc. In our feasibility study we have assumed a 
spherical shape. A sketch of a Laue lens is shown in Fig. 3. It has the appearance of a spherical corona, with 
inner and outer radii which depend on the energy bandwidth of the lens (see below). Assuming a paraboloidal 
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Figure 2. Peak reflectivity vs. energy for HOPG (002), Cu (111) and Ge (111). Tlie crystal thickness is 0.5 cm for the 
HOPG and 0.2 cm for both Copper and Germanium. A value of 70 = 1 and 77 = 1 arcmin is assumed. 




Figure 3. Left: Sketch of a Lauc lens with the shicklcd detector in its focus. Paths of paraxial photons are also shown. 
Right: Top view of the lens with a crystal arrangement according to an Archimedes spiral. 

shape, assuming a spherical geometry the image of a source at infinity is at half of the sphere radius. Spherical 
geometries are the most used also to focus diverging X-ray beams. ^ 

3.2. Crystal orientation 

Given that flat crystals arc more feasible, wc have assumed plane parallel crystal tiles of sizes , d± and t (see 
right panel of Fig. 4), where is the length of the crystal side parallel to the axis f of the mean lattice plane, 
d± is the length of the crystal side perpendicular to f , and t is the crystal thickness. The lenght of the crystal 
sides and d± should be as small as possible in order to better approximize the lens shape. In fact their sizes 
are a compromise between the focusing capabilities of the lens and number of crystal tiles required. 

Before discussing the arrangement of the single crystals in the lens we introduce a reference frame that 
analytically defines the position and orientation of each crystal (see Fig. 4). Given the spherical shape of the 
lens, a cylindrical reference frame, with the z axis coincident with the lens axis and its origin O in lens focus, 
it is natural choice. In this reference frame the position of a generic crystal is singled out by means of three 
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Figure 4. Left: Reference frame for the lens crystal position and orientation. Right: orientation of the lattice planes of 
each crystal tile. In order that paraxial photons of energy E which impinge on the crystal center C are reflected along 

the direction CO, the crystal must be rotated around the y axis by an angle 6b, where 6b is related to E through the 
Bragg law. The point C is also used to identify the position of the crystal tile. 



coordinates: the coordinates (vc, Zc) of the crystal center C, where Vc is the distance from the lens axis and Zc is 

the distance from the focal plane, and the rotation angle 9 of the crystal around the y axis (angle between the 
crystal axis and the z axis). To correctly focus photons of energy centroid E, it is needed that 9 = 9b (Bragg 
angle corresponding to E, see Fig. 4), and that 

Zc = F{2cos9b-1), (8) 

Tc = Zc tan 26*5 = 2F sin 6*5 . (9) 

where F is the focal distance of the lens (i.e., half radius of the spherical corona radius). From the previous 
equations, the rotation angle, and so the energy of the diffracted photons, depend on the crystal position through 
the relation Vc/zc = tan 20b. 

3.3. Lens energy passband and projected area 

The nominal energy passband [Emim Emax) of the lens is related to the values of the outer and inner radii Tmax 
and rmin-, respectively, of the spherical corona, through Eq. (9) and the Bragg law. We get 

Emin = -j (10) 

(*hkl ^max 

and 

_ hcF 

^hkl "^vnin 

One of the key parameters of a Laue lens is the projected area of the lens in the focal plane. If the entire 
spherical corona was entirely covered with crystals, the projected area would be that of circular corona, Acorona = 
'^(''^max " '''min)- ^^^^ area can be maximized in two ways: or increasing the lens focal length, or, for a given F, 
using materials with lower dhki- It can be easily shown that Acorona oc F"^. We notice however that higher focal 



lengths require smaller spreads of the mosaic crystals to avoid a larger dispersion of the diffracted photons in 
the focal plane (Point Spread Fimction, PSF). The PSF linearly increases with F. 

For a given F and lens energy passband, r increases with Ob, and thus, through Eq. (2), the lens projected area 
can be increased using materials with lower dhki- In principle dhki could be lowered even using highest diffraction 
orders, but this would imply lower diffraction efficiency, since Q{E) (see Eq. 6), through Fhki, decreases with the 

diffraction order. 

If the spherical corona cannot be completely covered with crystals (this is the case in which flat crystal tiles 
are available), it is useful define a filling factor f, which gives the ratio / = Acrys/^corona, where Acrys is the lens 
area covered with crystals. It can be easily shown that, for a given lens energy passband, / increases with the 
ratio between the focal length and the crystal tiles dimensions and thus using a greater number of tiles, which 
allows to obtain a better approximation of the spherical corona. 

3.4. Arrangement of the crystals in the lens 

While the lens energy passband, for a given F, depends on the inner and outer radii of the spherical corona, it 

can be shown that the lens effective area (defined as the projected area of the lens in the focal plane times the 
reflection efficency) depends on energy in a way which is related to the arrangement of the crystal tiles in the 
lens. 

In order to get an energy behavior of the lens effective area as smooth as possible, we followed the suggestion 
by Lund^ of arranging the crystal tiles according an Archimedes spiral: 

ri = ro + (d||+25)^, (12) 

where rj is the distance of the crystal i from the lens axis, dy is the length of the crystal tile along the lens radial 
direction (see Fig. 4), 5 is the extra-space needed for its accommodation, and is the azimuthal angle of the 
crystal center i, i.e. the angle between the radial vector fi of the crystal i and the radial vector ro of first crystal 
positioned in the lens. The coordinate Zc is set according to the already seen relation, Eq. (8). 

So, for a given F, the crystals are positioned with the following procedure: once the Bragg angle = ^min is 
established on the basis of the desired high energy threshold of the lens (see Eq. 11), the first crystal is positioned 

with its center at the coordinates (rg, zq), where tq = zq tan {29o) and zq = F(2cos9q — 1), and azimuthal angle 
(arbitrary position). Then, the coordinates of the successive crystal are determined by increasing the angular 
coordinate by Acf): 

A(Ai = - = 2arctan ( - — -^—r—^], (13) 

\2ri-i - d\\+d J 

that guarantees the needed room between contiguous crystal tiles. 

Once (pi is determined from Eq. 13, we can use Eq. (12) to obtain ri, and then the Bragg angle 9i and the 
distance of the crystal from the focal plane. 

Starting from the innermost crystal, all the other crystals are recursively positioned in the lens following the 
above procedure until the angle 6i of the last positioned crystal tile does not achieve a value 9max corresponding 
to the low nominal energy threshold of the lens Emin (see Eq. 10). 

3.5. Crystal positioning accuracy 

One of the main tasks to be faced in the integration of many (hundreds to thousands) mosaic crystals for a Laue 
lens is the accuracy in their orientation. Indeed this is crucial to reflect photons from the numerous crystal in 
the same focal plane position. The maximum error A9 in the nominal value of Oi shoiild be a fraction of the 
spread /? of the crystals. This angular error is due to the uncertainty in the knowledge of the crystal optical axis 
and to the tolerance in the machining of the lens support. An error A6 in the orientation of a crystal produces 
a linear displacement of the focal plane centroid of the crystal PSF by Ar = 2FA6. 



Table 1. Base configuration used in calculations. 



Mat(Tial 


Copper 


Miller Iiulicc^s 


(111) 


Energy Band 


60-200 keV 


Focal Length 


300 cm 


Tile Dimensions 


1.0 X 1.0 X 0.3 cm^ 


Mosaic Spread (FWHM) 


5 arcmin 


Number of Tiles 


2368 


Inner /Outer Radius 


8.92/29.7 cm 


Filling Factor 


0.94 



Table 2. Dependence of the tiles number, the weight of the crystals, the inner and outer radii {vmin and Vmax) and the 
filling factor for various focal lengths (F). 



F (cm) 


Tiles Number 


Weight (kg) 




rmax (cm) 


Filling Factor 


200 


1038 


2.78 


5.95 


19.8 


0.926 


300 


2368 


6.34 


8.92 


29.7 


0.939 


400 


4238 


11.3 


11.9 


39.7 


0.940 



4. DEVELOPED SOFTWARE 

For our feasibility study, we developed a software code (SWC) which analytically describes the lens geometry, 

its focal length, the mosaic crystal properties, the crystal sizes and their nominal orientation, the shape and 
size of the expected photon spot in the focal plane due to a paraxial beam of incident photons, the intensity of 
the diffracted beam produced by all crystals, the effect of misalignment of the mosaic crystals in the lens. All 
the lens parameters and crystal properties can be set at the time of run of the SWC or before according to our 
needs. For each set of given parameters, the SWC first computes the arrangement of the crystal tiles and then 
the expected effective area of the lens as a function of energy for on-axis incident photons. 

The initial development of the SWC was done using the proprietary programming language rsi-IDL, but we 
are porting it to the Python programming language. We decided to port the SWC because Python is an open 
source highly object oriented program language with a plenty of scientific libraries, and it is largely used by the 
scientific community. 

Usually the IDL software is used via a macro of commands which contains the complete set of instructions 
to do all the basic calculations and to give an initial output which is elaborated again through Python scripts. 

5. EXPECTED PERFORMANCE RESULTS 

We report here a few relevant results of our feasibility study. We have assumed, as a reference, a spherical lens 
with parameters summarized in Table 1. It makes use of fiat mosaic crystals of Cu (111) with a cross section 
of 10x10 mm^ and mean lattice plane perpendicular to the cross section. The nominal energy band of the lens 
is assumed to be 60 to 200 keV. We have computed the behavior of lens effective arcia Ag// with energy for 
different values of the focal length, the crystal thickness, the crystal mosaic spread and the uncertainty in the 
orientation of the crystals. The results are shown in Fig. 5. For the uncertainty in the orientation of the crystal 
we have assumed a Gaussian distribution with a FWHM of 12 arcmin. 

In Table 2 wc show the dependence of the crystal tiles number, lens crystal weight, inner and outer lens radii 
and filling factor as a function of the focal length, for three different values of F (200, 300, 400 cm). 

As can be seen from Fig. 5 (panel D), the effect of the crystal misorientation does not significantly influence 
the effective area. However it has a strong impact on the lens PSF and thus on the lens focusing factor given by 




Figure 5. Expected on-axis effective area with energy for different lens parameter values. Panel a: dependence on the 
focal length; panel b: dependence on the mosaic spread; panel c: dependence on crystal thickness; panel d: dependence 
on crystal orientation uncertainty (see text). 

G = eofph^eff /^d^ , where is the detection area which collects the fraction fph of reflected photons, and e/j 

is its efficiency. The focusing factor is a key parameter of the lens: its flux sensitivity depends on G. In Fig. 6 
we show the focusing factor for the reference lens configuration with parameters given in Table 1 for a nominal 
orientation of the crystal tiles in the lens, and for crystal orientations which are distributed around the nominal 
values according to a Gaussian fimction with a FWHM of 12 arcmin. The detector area Ad has been chosen in 
order to get the maximum G: in the first case Ad = 1.13 cm'^ with a collection fraction fph ranging from 0.60 to 
0.74, while in the latter case Ad = 1.77 cm^ with fph from 0.48 to 0.60. As can be seen, a sensitive reduction of 
the focusing factor is apparent when the crystals are not correctly oriented in the lens. 

6. PROSPECTS 

The results obtained by our feasibility study are very promising for the prospects of the hard X-ray astronomical 
observations (> 70 kcV). The Lauc lenses appear to be the most promising focusing optics to overcome the 
sensitivity limitations of the current hard X-ray instruments. A key advantage of the Laue lenses is their very 
low weight (see Table 2), much lower than the multilayer optics now under development for efficiently focusing 
photons with energies below 70 keV. 

The next step of our project is the development of a lens demonstration model in which we get ready the lens 
crystal assembling technique. Soon after we intend to develop a prototype model to be tested also in a balloon 
flight. 

We have proposed to use Laue lenses for a high energy experiment aboard a long duration balloon flight at 3 
mbar of residual atmospheric pressure. However the best exploitation of the Laue lenses is in telescopes aboard 
satellite missions in which high focal lengths are possible. The use of satellite formations, with a satellite with 
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Figure 6. Focusing factor for an ideal lens, and for a lens with flat mosaic crystal tiles misoriented with respect to the 
nominal rotation angle according to a Gaussian function with a FWHM of 12 arcmin (see text). The detection efficiency 
is assumed to be to = 1- 

aboard the lenses and a slave satellite at a distance of tens of meters with aboard the focal plane detectors should 
be ideal for exploiting the capabilities of the Laue lenses. Such satellite formations are included in the strategic 
program of ESA. 
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